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Abstract 

The mapping of a rotational dynamics on a harmonic oscillator one is considered. The 
method is used for studying the stabilization of the rigid top rotation around the intermediate 
moment of inertia axis by orbiting particle. 


1 Introduction 

The quantum rotation is a specific type of excitation of microscopic system: hadrons, nuclei, 
molecules, and even atoms. The rotational excitations of molecules and nuclei have been studied 
in more detail. Electronic excitations are much higher than vibrational ones for most so-called 
” normal” molecules. Therefore they may be described adequately in the Born-Oppenheimer ap 
proximation. There is no analog of the Born-Oppenheimer approximation for atomic nuclei. Yet 
the occurrence of the rotational bands with strong (nearly 100 single particles) E2-transitions 
between neighboring states shows the existence of the collective rotation. All nucleons participate 
cooperatively into this collective motion with internal degrees of freedom being frozen completely 
or partly. The rotational excitations are grouped into rotational bands having states characterized 
in simplest case of a rigid axially-symmetric top by the energy (h = 1) 

E = I(I+ l)/23, (1) 

and quantum number I of the total angular momentum. 3 is a moment of inertia. A simplest 
non-axial system is a rigid asymmetric top with the Hamiltonian 

H = A 1 I 2 l +A 2 I 2 2 + A 3 ll (2) 

where I a are the projections of the total angular momentum operator on the BFF (body-fixed 
frame) axes a — 1,2,3. The rotational constants A a = l/(23„) depend on the principal moments 
of inertia $ s a . The rotational band of an asymmetric top consists of rotational multiplets, i.e., of 
the levels with the same value of the quantum number I. Besides /, these levels are characterized 
by the irreducible representations of the group D 2 = {l, 3 ?i, 3 ff 2 1 ^ 3)5 which contains the identity 
operator and three 180° rotations around the BFF axes $ft a (7r) = exp(—inl a ). The irreducible 
representations of D 2 are labeled ai,a 2 ,bi,b 2 . They correspond to the eigenvalues r a = ±1 of the 
operators 9i Q . The subscripts 1 and 2 label even and odd symmetry levels with respect to the 3?2 
rotation, a and b label even and odd symmetry levels with respect to the 9ft 3 rotation. Note that 
riT 2 r 3 = 1 is true for each four representations. 
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Different methods are used to solve the Hamiltonian (2). As the first step we will calculate 
the classical rotation energy E as the function defined in the system phase space (a rotational 
energy surface). The phase space of a rotational motion is formed by three Euler angles <J>, 
and three conjugated momenta p<t,,p#,P 4 ,- The absolute value of the angular momentum I and its 
projection I z = p<j, on the z-axis of the space-fixed frame are the integrals of motion. It is suitable 
to do the canonical transformation [ 1 ] to new conjugated variables I and qi , I z and q z , I 3 and q. 
Since qi and q z are cyclic variables, the phase space of a rotational motion is two-dimensional. It 
is convenient to map it on the surface of the sphere of the / radius with a center in the origin 
of the BFF. The point on the sphere with coordinates 0 and y? determines the orientation of the 
vector I in the BFF. The canonical transformation enables us to relate the conjugated variables I 3 
and q to the angles 6 and y?. For I z = qi = 0 and an arbitrary q, we have cos 6 = / 3 / / , ip = f - q. 
Thus, the trajectories of the tip of I on the phase sphere are classical trajectories of the system 
in its rotational phase space. When the rotational energy is close to A\I 2 or A 3 / 2 , where A\ and 
A 3 correspond to the smallest or the largest moment of inertia, the classical trajectories are small 
ellipses around axes 1 or 3. They represent precession motion around these axes. The trajectories 
close to axis 2 with the intermediate moment of inertia are unclosed. They represent unstable 
motion. A small deviation from this axis takes a top away from it. 


2 Precession motion 

Let us begin with classical precession. It is described by the Euler equations 

i a = {H , /*}, a = 1,2,3, (3) 

for the projections of the vector I on the BFF axes. In this equation {...} are the Poisson 
brackets. Let us introduce the classical concept of stationary rotation axis defined by the three 
equations {H, I 0a } = 0. The stationary states I 0 are identical with the fixed points of the energy 
surface. There are three stationary axes coinciding with the principal axes of a rigid top. For 
small precession around axis 1 (7) ~ /; I 2 , 1 3 <C I) Eqs. (3) have the form 

I 2 = ~{A 3 - Ai)II 3 , 

h = (A 2 - Ax)// 2 . ' 

They describe a harmonic oscillation motion 

/ 2 (f) = loyA^/Ai — 1 coswnt, Iz(i) = —iQ\jA 2 l A\ — lsino>nt, (5) 

with small amplitude io and frequency 

u>n = 2IyJJ~A 2 — Ax)(A 3 — ~A\). (6) 

In the BFF, the I vector precess around axis 1 and, in the laboratory frame, the top precess around 
the angular momentum I. It follows from the stability condition u 2 > 0 that the precession motion 
around the axis with the smallest or the largest moment of inertia is stable and around the axis 
with the intermediate one is unstable. 

To obtain the energy level structure of a rotational band corresponding to a precession motion 
one must quantize this motion. It can be done by different methods. 


130 



Holstein- Primakoff Representation of Angular Momentum Operators. The method developed 
by Marshalek for a quantized rotator [2] is based on the one boson realization of angular momentum 
operators [3] 

7 1 = 7-6+6, 7+ = h + ih = &V2/ - 6+6 = (/_)+, (7) 


in the space of wavefunctions 


s piv — 2 ^ 
K=-l 


[(7 — Ky.yi 2 


o), 


( 8 ) 


where v is the quantum number of a state in a rotational multiplet. The state | 0) corresponding 
to K = 7 is a vacuum state of the boson creation and annihilation operators 6+, 6. For describing 
precession motion we expand the square-root operator in 7+ and 7_ in the series of a small quantity 
n/7, where h = 6+6 is the boson number operator. In the harmonic approximation (n < 7), the 
Hamiltonian (2) is quadratic in boson operators 

H = AJ 2 + A 2 + A 3 ) + I(A 2 + A 3 - 2A 1 )b + b+ ! -{A 2 + A 3 )(b+b + 66). (9) 

The Hamiltonian (9) can be diagonalized by a canonical transformation 

6 = u/3 -f u/?+, | u | 2 — | v | 2 = 1, (10) 

to new boson operators /3 and /?+ . The energy of the lowest levels is given by 

Ej n = AiI(I A 1) +uqi(n + 1/2), n = 0,l,2, ... . (H) 


The quantum 
function 


number n describes the precession motion. For the state with n — 0 
Y?/o =| u | _1/2 exp (^ 6+6+ ) I °)> 


the wave 

(12) 


localizes near rotational axis 1. It corresponds to the sharply localized orientation of the angular 
momentum I along the positive direction of axis 1. Eq. (12) is not the eigenfunction of D 2 
symmetries. Thus, the harmonic approximation is a ” symmetry-breaking approximation.” Being 
a linear theory it cannot describe tunneling the vector I through a potential barrier separating 
two degenerate minima. 


Bargmann Representation [4]. Let us consider the complex variable 


c 


= x + iy = e t( ^ 



(13) 


which is the stereographic projection of a point on the pha.se sphere with polar angles 0 and y? 
on the plane passing through the south pole. As shown in Ref. [5], it is possible to construct the 

SU{ 2) group representation in the space of the polynomials 


v>(C) = E «a-c /+ k '- 

K——I 


(14) 
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The angular momentum operators in this representation have the form 

h = I( + “ C 2 )^ I 2 = + 2^ + + C^- (15) 

Many problems in nuclear and molecular physics can be treated by this approach [6]. With the 
operators (15), the Schrodinger equation for a top is reduced to the Heine equation [7]. To describe 
for example the precession motion around axis 3 we should consider the approximation | ( |<C 1. 

Approximate Solution of Reccurence Relations. This method is based on the approximation of 
recurrence relations by a second-order differential equation for high-I values [8]. The eigenfunction 
of the Hamiltonian (2) can be written in the form 

/ 

^ IMv = a iKvD l MK {'d ), (16) 

K=—[ 


where M and K are the quantum numbers of the operators I z and / 3 respectively. D I MK is 
the Wigner function depending on the Euler angles d. For the coefficients a/^v, the three-term 
recurrence relation is obtained 

Pk'ai.K- 2,u + (Hkk — Eiv)aiKv + Pk+2<ii,k+2,is = 0 , ( 17 ) 

where 

Hk = Hkk = \{A\ + A 2 ) + |(2A 3 — Ai — A 2 )K 2 , 

Pk = Hkk-2 = \(A X - A 2 )[(I - K + 2 )(/ + K - 1)(/ - K + 1 )(/ + A')] ,/2 . ( 18J 

By using the small parameter 6 = [1(1 + l)] -1 / 2 let us introduce the continuous variable k — KS. 
We will treat the coefficients P and H as the smooth functions of this variable. As a result, the 
recurrence relation (17) may be rewritten in the form of the Schrodinger equation 


/(/ -(- l)Ha/„(fc) = Ei v ai v (k ), 


(19) 


with the Hamiltonian 


H = [P(k + 26) + P(k)] cos 2 pS + i[P(k A 26) — P(k )] sin 2 p6 + H(k), (20) 


where p = —id/dk is the canonically conjugated momentum to the coordinate k. 

In the harmonic approximation K <C /, Eq. (19) is reduced to the second-order differential 
equation 

^ + 2m[E - AJ(I + 1) - V(K)]a Iu = 0, (21) 

in the space of quantum numbers K . The equation describes the motion of the angular momentum 
I with effective mass m = 1/[2(A 2 — Ai)/(/ + 1)] in potential V(K) = (A 3 — Ai)K 2 . This harmonic 
oscillator has the energy spectrum (11) and the wave function 


a>in — 


2 

1 2 w 

2 n n! \ 

' 7r(2/+ 1) _ 


H n 



(22) 
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where w = [{A 3 - Ai)/(A 2 - ^i)] 1/2 , and H n is the Hermite polynomial. We have again a sharp 
localized state with broken symmetry. 

The rotational dynamics of real many-body systems is more complicated than that of a rigid 
top because of the centrifugal and Coriolis forces. However the harmonic approximation can 
be used in this chase to understand the physical nature of the phenomenon under study. The 
example is bifurcations in quantum rotational spectra considered in Ref. [9]. In the next section 
we consider another problem having a bearing on a precessional motion. 


3 Precession Motion around Intermediate Moment of In- 
ertia Axis 


As we proved above, the rotation of a rigid top around the intermediate axis is unstable. This is 
not a case if a system has additional degrees of freedom apart from rotational ones. We are going 
to consider a particle coupled with an asymmetric rigid top. This is the situation of one-electron 
Rydberg states in triatomic molecules, such as H 2 0 [10], and rotational bands in odd deformed 
nuclei [ 11 ]. 

The Born-Oppenheimer approximation breaks down in the molecular Rydberg spectra. As the 
total angular momentum / increases, the rotational levels pass from Hund’s case (a) or ( b ) (the 
strong-coupling scheme), where the electronic splitting is large compared with the rotational one, 
to the Hund’s case (d) (the weak-coupling scheme), where it is small. The model of an isolated j- 
complex is widely used for the description of the transition from strong to weak coupling in nuclear 
[ 11 ] and molecular [12] physics. This approximation is valid if coupling the states of a j-complex 
with other Rydberg states is small compared to the Coriolis coupling. The assumption means 
that one-particle angular momentum j is an integral of motion. As / increases, the momentum 
j decouples from the molecular ion core and couples to the axis of rotation with the maximal 
moment of inertia. 

Let us consider the effective Hamiltonian describing the two degrees of freedom: rotational and 
one-particle. The rotational part of the Hamiltonian is the kinetic energy of a rigid top. We will 
use the self-consistent field approximation for describing one-particle motion. The non-sphencal 
part of this field can be written in terms of the particle multipole momenta <?v( r) as follows: 
y — Hxti <?A^<?A M (r), where A-values are even for the reflection symmetric field. For an isolated 
j-complex, the one-particle part' of the Hamiltonian can be expressed in terms of the spherical 
tensor operators 7\„(j). Thus, the effective Hamiltonian of the system for an isolated j-complex 
is algebraic with symmetry SU{ 2) ® SU( 2). We consider this Hamiltonian in the quadrupole 
approximation when the components of mean field V with A > 2 are smaller than the ones with 
A = 2. It is convenient to write the Hamiltonian in the coordinate system fixed by the principal 


inertia axes: 

H = £[A q (/ 0 - j a ) 2 + 9ajl\, 

a 


(23) 


where gi = - g 2 = {go sin 7)/\/3, g 3 = g 0 cos 7 (-00 < g 0 < + 00 , 0 ° < 7 < 60°) are the parameters 
of the quadroupole part of the self-consistent field, which is diagonal in the considered frame. In 
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classical mechanics, the system dynamics is described by the equations of motion: 

ja 2 €otj3~f[ — Ap^Ip jp) -h ff0^0\jy ? 


( 24 ) 


where repeated indexes are summed. 

The stationary state Io and jo of Eqs. (24) is determined by the eight algebraic equations 
since the two integrals of motion I 2 and j 2 exist. Therefore, the stationary state with an arbitrary 
orientation does not exist. There are two types of stationary states of Eqs. (24), corresponding to 
the lowest level of the rotational multiplets. The three aligned states S Q : I a = /, j a = j, a = 1, 2, 3 
with the energy 

E a = A a (I - j) 2 + g a j 2 , (25) 

have the parallel vectors Io and j 0 aligned along axis a. In the three plane stationary states S a p, 
these vectors are placed in the (a/?)-plane. As 7 increases the sequence of stationary states S a and 
S a p with the minimal energy for given 7 leads to the aligned state with the maximal moment of 
inertia. The transition from aligned state S a to plane one S a p is accompanied by the bifurcation 
of the C^-type [9] at the critical angular momentum 


7 a0 — J 


1 + 


90 


2A a 


+ 


( 9a ~ 90 \ 
V 24. ) 


2A 


(9c, — 9p)Ap 

A a (A a - Ap)\ V 


(26) 


The index a denotes the axis from which the angular momenta Io and j 0 decouple, while index 
/? denotes the axis to which they approach. Both indexes (a/3) denote the plane, in which these 
vectors move for j < I a p < I < Ip a . Another pair of critical momenta determines the similar 
bifurcation points in the region I < j. The expression for these values is the same as Eq. (26) 
except the sign before a square root. We will consider below only the region 7 > j in the case 
when Ai<A 2 < A 3 , g 0 > 0. 

The precessional motion near the stable stationary state is described by linearized set of Eqs. 
(24). Four linear differential equations describe two normal modes corresponding to the small 
harmonic oscillations of vectors I and j. The frequencies of these modes are obtained as the roots 
of a biquadratic equation. We begin our consideration with the precession near the 5'1-state in 
the weak-coupling limit A a I » g a j . To simplify expressions let us use assumption I»j. The 
precessional mode with smaller frequency u> u (6) represents the precession of the vectors I and 
j with different amplitudes. While the I vector circumscribes according to Eqs. (5) an elliptical 
cone around axis 1 with the amplitude «o, the amplitude of the j vector is j/I times smaller that 
of the I one. Thus, this mode involves the precession of core angular momentum vector R = I — j, 
which coincides, in considered approximation, to rigid top precession. Another normal mode with 
the frequency = 27 A\, which is equal to the core angular velocity around axis 1, involves the 
uniform rotation of the angular momentum vectors around this axis with identical amplitudes: 


h(t) =h(t) = *o cos u u t, h(t)=j 3 (t) = *o sin uji 2 t. (27) 

Consequently the momentum R does not participate in this motion. 

Now we consider precession near the stationary state S 2 , corresponding to the central axis. An 
orbiting particle stabilizes this state due to the anisotropic (quadrupole) interaction with a top. 
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The small values of ( A 2 — Ai)/A 2 favors stability of the aligned state S 2 lowering its energy E 2 
relative to E\ (see Eq.(25)). To clarify the physical picture we will use the approximations I»j 
and (A 2 — Ai)/A\ « 1 in describing precessional motion. The smaller frequency 


U> 2 i 


2 


(A 2 — y4i)(A3 — A 2 )(I 2 1 



(28) 


vanishes at the critical momentum I 2 \ (26). The time-dependence of the angular momentum 
components for this mode is defined by 

I\ (t) = i 0 [A 3 /A 2 - l] 1/2 cos w 21 t, j29) 

I 3 (t) = -i 0 [(A 2 /A 1 -\)(I 2 JP - l)] 1/2 sin u 2l t, 

and j(i) = (j/I)l(t). The interpretation of this result is straightforward. Since the amplitude of the 
j-vector is small compared to that of the I ones, the considered mode represents the precession of 
the core angular momentum R similar to the rigid top precession (5). Yet unlike the latter, the tips 
of vectors I and j move on elliptic orbits stretched along axis 1 if the angular momentum I is close 
to I 2 i. This is just a consequence of the bifurcation, which shifts the angular momentum vectors 
into the (12)-plane. For another normal mode of frequency u 22 = 2 M 2 , the time-dependence of 
angular momentum components has the same form as in Eq. (27). 

Thus, we have shown that the precession motion around the axis with intermediate moment of 
inertia is possible for a system consisting of a particle anisotropically coupled with an asymmetric 
top. The isolated ^-complex approximation is used in considering this phenomenon, lo examine 
it one can solve the classical equation of motion for a particle coupled with an asymmetric top 
without this restriction. The equations involve two different time scales, fast particle motion and 
slow core rotational motion. After averaging on the fast motion, one can obtain the closed set of 
nine equations for components of I and particle angular and quadrupole momenta, d he averaged 
equations can describe the stabilization phenomenon and the precession around intermediate axis. 
This insight into the problem reveals the close relationship of our stabilization effect with that of 
the Kapitza pendulum [13]. Another intriguing analogy is the new discovery in planetary science 
where it has been shown recently [14] that the Moon stabilizes the chaotic wobble of the Earth s 
rotational axis, which is unstable due to orbital coupling with other planets. Thus, without the 
Moon, large variations in obliquity resulting from the chaotic wobble might have driven dramatic 
changes in the Earth’s climate. There are two fundamental distinctions of our problem from 
considered above: we deal with the isolated and quantum system. 

The above found precessional frequencies are associated with the splitting between the lowest 
levels of a multiplet. To obtain this result one must quantize the precession motion. It can 
be done by using the Holstein-Primakoff representation. The result for lowest multiplet levels 
corresponding to the stationary state S a is 

E In = E a (I) + w a i(m + 1/2) + u a2 {n 2 + 1/2) (30) 

where E a is given by Eq. (25), and n x and n 2 are the numbers of bosons in corresponding mode. 
The boson operators b% and b k connect the odd and even with respect to 'the (^-transformation 
states inside a rotational multiplet. Consider, for example, the precession around axis 1 in the 
weak-coupling limit. In this approximation, any rotational multiplet (/-multiplet) consists of R- 
multiplets with the quantum numbers R = I — j + m, m = 0, 1, ...,2 j. The levels in /Cmultiplets 
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with the same quantum number R but different I are degenerated. The frequency (6) describing 
the precession of the R-vector is equal to the splitting between the lowest levels ai and b\ (or 
a 2 and b 2 ) in a R-multiplet. Another frequency uq 2 = 21 A\ is equal to the splitting between 
the lowest levels of adjacent /2-multiplets belonging to the same /-multiplet. The situation for 
the precession around the intermediate axis is more complicated. According to the precessional 
approximation, the lowest states of a multiplet involve two groups of roughly equidistant levels, 
which are described by Eq. (30). But a smaller frequency vanishes in a critical point and the 
precessional approximation becomes inappropriate in this region. This means the redistribution of 
multiplet levels, which provides a method for the identification of the intermediate axis precession 
in an experiment. 
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